Crossover between two different magnetization reversal modes 
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The magnetization reversal in ordered arrays of iron oxide nanotubes of 50 nm outer diameter 
grown by atomic layer deposition is investigated theoretically as a function of the tube wall thickness, 
d w . In thin tubes (d w < 13 nm) the reversal of magnetization is achieved by the propagation of a 
vortex domain boundary, while in thick tubes (d w > 13 nm) the reversal is driven by the propagation 
of a transverse domain boundary. Magnetostatic interactions between the tubes are responsible for 
a decrease of the coercive field in the array. Our calculations are in agreement with recently reported 
experimental results. We predict that the crossover between the vortex and transverse modes of 
magnetization reversal is a general phenomenon on the length scale considered. 

PACS numbers: 75.75.+a, 75.10.-b 



I. INTRODUCTION 

Magnetic nanoparticles, particles of nanometer size 
made from magnetic materials, have attracted increasing 
interest among researchers of various fields due to their 
promising applications in hard disk drives, magnetic ran- 
dom access memory, and other spintronic devices Be- 
sides, these magnetic nanoparticles can be used for poten- 
tial biomedical applications, such as magnetic resonance 
imaging (the nanoparticles can be used to trace bioana- 
lytes in the body), cell and DNA separation, and drug 
delivery!^ To apply nanoparticles in various potential de- 
vices and architectures, it is very important to control 
the size and shape and to keep the thermal and chemical 
stability of the nanoparticles!^ 

The properties of virtually all magnetic materials are 
controlled by domains - extended regions where the spins 
of individual electrons are tightly locked together and 
point in the same direction. Where two domains meet, a 
domain wall forms. Measurements on elongated magnetic 
nanostructures^ highlighted the importance of nucleation 
and propagation of a magnetic boundary, or domain wall, 
between opposing magnetic domains in the magnetiza- 
tion reversal process. Domain-wall propagation in con- 
fined structures is of basic interest P^I p or instance, by 
equating the direction of a domain's magnetization with 
a binary or 1 , a domain wall also becomes a mobile edge 
between data bits: the pseudo-one-dimensional structure 
can thus be thought of as a physical means of transport- 
ing information in magnetic form. This is an appealing 
development, because computers currently record infor- 
mation onto their hard disk in magnetic formP^ 

The trusty sphere remains the preferred shape for 
nanoparticles but this geometry leaves only one surface 
for modification, complicating the generation of mul- 
tifunctional particles. Thus, a technology that could 
modify differentially the inner and outer surface would 
be highly desirable.^ On the one hand, over the past 
years there has been a surge in research on nanocrystals 



with core/shell architectures. Although extensive stud- 
ies have been conducted on the preparation of core / shell- 
structured nanoparticles, the fabrication and characteri- 
zation of bimetallic core/shell particles with a total size of 
less than 10 nm and with a monolayer metal shell remain 
challenging tasksP^ On the other hand, since the discov- 
ery of carbon nanotubes by Iijima in 1991 pi intense at- 
tention has been paid to hollow tubular nanostructures 
because of their particular significance for prospective 
applications. In 2002 Mitchell et a/P^ used silica nan- 
otubes offering two easy-to- modify surfaces. More re- 
cently, magnetic nanotubes have been growrl^^ that 
may be suitable for applications in biotechnology, where 
magnetic nanostructures with low density, which can 
float in solutions, become much more useful for in vivo 
applications.^ In this way tiny magnetic tubes could pro- 
vide an unconventional solution to several research prob- 
lems, and a useful vehicle for imaging and drug delivery 
applications. 

Although the magnetic behavior of nanowires has been 
intensely investigated, tubes have received less attention, 
in spite of the additional degree of freedom they present; 
not only can the length, L, and radius, R, be varied, 
but also the thickness of the wall, d w . Changes in thick- 
ness are expected to strongly affect the mechanism of 
magnetization reversal, and thereby, the overall magnetic 
behavior! 20 * 21 * However, systematic experimental studies 
on this aspect were lacking for a long time, mostly due 
to the difficulty in preparing ordered nanotube samples 
of very well-defined and tunable geometric parameters. 

We recently reported the synthesis and magnetic char- 
acterization of a series of Fe304 nanotube arrays (length 
L = 3 /im, radius R = 25 nm, center-to-center distance 
D = 105 nm, and wall thickness 2.5 nm < d w < 22 
nm), prepared by atomic layer deposition (ALD) in a 
porous alumina matrixP^In this series, the magnetic re- 
sponse of the array, characterized by the coercive field 
H c and the relative remanencep^ vary strongly and non- 
monotonically as a function of d w . For the thinner tubes, 
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H c is enhanced by increasing d wi until d w — 13 nm, 
where it presents a maximum of about 780 Oe. For fur- 
ther increases of d wi the coercive field decreases. A quan- 
titatively similar behavior was also observed in Ni8oFe2o 
nanowire arrays, 24 a different system in terms of geome- 
try, material and preparation techniques. 

This convergence of experimental observations may 
reflect an underlying general phenomenon. Therefore, 
this paper focuses on the investigation of the non- 
monotonic behavior of the coercive field in ferromag- 
netic nanotube arrays, a question that has remained un- 
explained until now. We start by modelling the mag- 
netization reversal and calculate H c for the system re- 
ported experimentally,^ then generalize our conclusions 
and quantitatively predict trends for other geometries 
and materials. 



II. EXPERIMENTAL METHODS 

Our approach to the preparation of magnetic nan- 
otubes of well-controlled and tunable geometric parame- 
ters and arranged in hexagonally ordered, parallel arrays 
is based on the combination of two complementary as- 
pects, namely (i) the utilization of self-ordered anodic 
alumina (AA) as a porous template, and (ii) the confor- 
mal coating of its cylindrical pores with thin oxide films 
by atomic layer deposition (ALD). 

Anodic alumina is obtained from the electrochemical 
oxidation of aluminum metal under high voltage (usu- 
ally 20 to 200 V) in aqueous acidic solutions! 25 1 26 1 Un- 
der certain proper sets of experimental conditions (na- 
ture and concentration of the acid, temperature and ap- 
plied voltage), the electrochemically generated layer of 
alumina displays a self-ordered porous structure. Cylin- 
drical pores of homogeneous diameter are thus obtained, 
with their long axis perpendicular to the plane of the 
alumina layer and ordered in a close-packed hexagonal 
arrangement. With our method, anodization of Al in 0.3 
M oxalic acid under 40 V at 8°C yields pores of ~ 50 
nm outer diameter and with a center-to-center distance 
of ~ 105 nm (an approach which we will call Method 
A); anodization in 1% phosphoric acid under 195 V at 
0°C yields pores of ~ 160 nm outer diameter and with a 
center-to-center distance of ~ 460 nm (Method B). 

Atomic layer deposition is a self-limited gas-solid 
chemical reaction.^ Two thermally stable gaseous pre- 
cursors are pulsed alternatively into the reaction cham- 
ber, whereby direct contact of both precursors in the 
gas phase is prevented. Because each precursor specif- 
ically reacts with chemical functional groups present on 
the surface of the substrate (as opposed to non-specific 
thermal decomposition), one monolayer of precursor ad- 
sorbs onto the surface during each pulse despite an ex- 
cess of it in the gas phase. This peculiarity of ALD 
makes it suitable for coating substrates of complex ge- 
ometry (in particular, highly porous ones) conformally 
and with outstanding thickness control.^ We have suc- 



cessfully used ALD to create Fe203 nanotubes in porous 
anodic alumina templates from two different chemical re- 
actions with similar results. In Method /, oxidation of fer- 
rocene (also called bis (cyclopentadienyl) iron, usually ab- 
breviated Cp2Fe) with an ozone/dioxygen (O3 / O2) mix- 
ture at 200° C yields a growth rate of ~ 0.2 A per cycle. 
Method II consists in the reaction of the dimeric iron (III) 
tert-but oxide, Fe2(O t Bu)6, with water at 140°C, with 
~ 0.25 A deposited per cycle. 22 Both methods yield a wall 
thickness distribution within each sample below 10%. 

The results are shown in Figs. 1 and 2. Smooth tubes 
of 50 or 160 nm outer diameters can be obtained, with 
aspect ratios on the order of 100. The thickness of the 
wall can be accurately controlled between 1 and 50 nm. 
Subsequent reduction of the Fe203 material by H2 at 
400° C results in the formation of the strongly magnetic 
phase Fes04, a transformation verified by X-ray photo- 
electron spectroscopy and accompanied by the expected 
color change from yellow, orange or brown (depending 
on the thickness) to black. The structural quality of the 
tubes is unaffected by reduction,^ a consequence of the 
very small volume contraction caused by it. Our ap- 
proach allowed us to systematically investigate the influ- 
ence of structure on magnetism in a series of samples of 
Fe304 nanotube arrays prepared according to Methods A 
and 77 and in which the wall thickness, d w , varies while 
all other geometric parameters are maintained constant. 
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FIG. 1. (Color online) Aspects of Fe2 03 nanotubes grown by 
atomic layer deposition (ALD) in a porous anodic alumina 
matrix, (a) Macroscopic view of the samples, consisting of 
a circular porous anodic alumina (AA) membrane containing 
the embedded tubes, surrounded by an outer circle of Al metal 
of 2 cm outer diameter; the tubes in the samples from left to 
right have walls of increasing thickness, approximately 0, 1, 
2, 4, 8, 12, and 16 nm, respectively (AA obtained according 
to Method A, ALD performed by Method I), (b, c) Scanning 
electron micrographs (SEM) of tubes embedded in the porous 
alumina matrix, observed in top view and at an angle at a 
break in the sample, respectively; the scale bar represents 200 
nm (AA obtained according to Method B, ALD performed by 
Method II). 

In a series of Fe304 nanotube arrays of varying wall 
thickness d w (all other geometric parameters being kept 
constant), investigated by SQUID magnetometry, we ob- 
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FIG. 2. Transmission electron micrographs (TEM) of isolated 
tubes from the alumina matrix, with walls of increasing thick- 
nesses, approximately of 1(a), 5 (b), and 13 (c) nm, respec- 
tively; the scale bar represents 50 nm (AA obtained according 
to Method A, ALD performed by Method II). 



served a significant dependence of the coercivity and re- 
manence upon the geometry. In particular, the coercive 
field H c can be tuned between and 800 Oe (0 and 80 
mT) approximately by properly adjusting d w . Most cu- 
riously, the dependence of H c on d w is not monotonic - 
H c reaches its maximum at d w ~ 13 nm and then de- 
creases for further increases of the wall thickness (Figure 
6b). We interpret this observation as arising from the 
coexistence of two distinct magnetization reversal modes 
in our system. Which of the two prevails in a given sam- 
ple is uniquely determined by the geometric parameters 
of the tube array. Thus, the cusp in the H c (d w ) curve 
corresponds to the crossover between the two modes of 
magnetization reversal. The following paragraphs detail 
the theoretical treatment of the two modes. 



III. TWO MAGNETIZATION REVERSAL 
MODES 

For isolated magnetic nanotubes, the magnetization 
reversal, that is, the change of the magnetization from 
one of its energy minima (M = Mo z) to the other 
(M = — M z), can occur by one of only two idealized 
mechanisms, the Vortex mode (V), whereby spins in ro- 
tation remain tangent to the tube wall, or the Transverse 
mode (T) , in which a net magnetization component in the 
(x, y) plane appears)^ In both cases, a domain boundary 
appears at one end of the tube and propagates towards 
the other, as illustrated in Fig. 3. Starting from the equa- 
tions presented by Landeros et a/p^ we can calculate the 
zero-field energy barrier as well as the width of the do- 
main boundary for each reversal mode as a function of 
the tube thickness, d w . Figure 4a, 4b, and 4c present our 
results for Fes04 nanotubes using M = 4.8 x 10 5 A/m 3 
and the stiffness constant A = 10 -11 J/mP2 Figure 4a 
witnesses a crossover at d w « 20 nm, showing that the 
V mode is more stable for thinner tubes, whereas thicker 
tube walls favor the T mode. This result can be qual- 
itatively explained as follows. A very thin tube should 
behave as a (rolled-up) thin film, in which the magnetic 
moments always tend to remain within the plane of the 
film. Conversely, tubes of large wall thicknesses approach 
the case of wires: surface effects are less crucial, but in- 
teractions between diametrically opposed regions become 
more important. 

The presence of a crossover in Fig. 4a allows us to 
expect a transition from the V to the T reversal mode 
with increasing values of d w . However, the curves can- 
not give the coercive field values directly because energies 
represent the difference between a completely saturated 
state and one with a domain boundary in the middle of 
the tube. Magnetization reversal, however, is initiated 
with a domain boundary at one end of the tube, a con- 
figuration that corresponds to a lower magnetic energy. 
From Fig. 4b we observe changes in boundary widths 
between 50 and 120 nm as a function of d w . A crossover 
is found, corresponding to the one that appears between 
the energy curves. We shall now proceed to calculate the 
switching field of an isolated magnetic nanotube assum- 
ing that the magnetization reversal is driven by means of 
one of the two previously presented modest 



A. Coercive fields 

For the T mode, the coercive field, if J, can be ap- 
proximated by an adapted Stoner-Wohlfarth modeP^ in 
which the length of the coherent rotation is replaced by 
the width of the domain boundary, wt (see Fig. 4b). 
Following this approach, 
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FIG. 3. (Color Online) Geometric parameters and magne- 
tization reversal modes in nanotubes. Arrows represent the 
orientation of magnetic moments within the tube. Left: A 
magnetically saturated tube (M = Mo z), with its geometric 
parameters, length L, radius R, wall thickness d w . Center: A 
tube during the magnetization reversal from Mo z to —Mo z 
via a vortex mode, V; the domain boundary, of thickness wv, 
migrates upwards. Right: The equivalent situation by means 
of a transverse reversal mode, T, with a domain boundary of 
thickness wt- 



where K (I) = |/i M 2 (1 - 3N Z (I)) and N z (I) corre- 
sponds to the demagnetizing factor along z, given by 
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For the V mode we use an expression for the nucle- 
ation field obtained by Chang et a/P^ When an exter- 
nal field with magnitude equal to the nucleation field is 
applied opposite to the magnetization of the tube, in- 
finitesimal deviations from the initially saturated state 
along the tube axis appear. The form of these devia- 
tions is determined by the solution of a linearized Brown's 
equation!^ Furthermore, it has been shown numerically 
that the solution of this Brown's equation is not a stable 
solution to the full nonlinear equation at applied fields 
larger than the nucleation field, and then the only possi- 
ble stab le sta tes are those with uniform alignment along 
the axis! 35 * 36 J Thus, the magnetization is assumed to re- 
verse completely at the nucleation field. For an infinite 
tube, the nucleation field for the V mode, H% , is given 

by 
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FIG. 4. (Color online) Energy barrier (a) and domain bound- 
ary width f&j as a function of tube wall thickness for the V 
(green) and T (purple) modes of magnetization reversal, (c) 
Coercive fields obtained from Eqs. (1) and (2) for the V and 
T modes, (respectively shown in green and purple) for tube 
wall thicknesses varying between 0.5 and 24.5 nm. The dashed 
lines of the curves have no physical meaning. We have used 
R = 25 nm 



Here J p (z) and Y p (z) are Bessel functions of the first and 
second kind, respectively. Equation (3) has an infinite 
number of solutions, and the physically correct solution 
is the smallest one. 

Figure 4c illustrates the coercive field of an isolated 
tube with d w varying from 0.5 to 24.5 nm. We can ob- 
serve a crossing of the two curves at d w = 13 nm ap- 
proximately, corresponding to a magnetization reversal 
for which both V and T mechanisms are possible at the 
same coercive field. At other given values of d Wl the 
system will reverse its magnetization by whichever mode 
opens an energetically accessible route first, that is, by 
the mode that offers the lowest coercivity. Therefore, the 
curve of coercivity vs. d w is the solid one, and the dotted 



5 



sections of curves have no physical meaning. Thus, our 
Fe304 tubes will reverse their magnetization by the V 
mode for d w < 13 nm, and by the T mode for d w > 13 
nm. Our calculations for an isolated tube reproduce the 
non-monotonic behavior of the coercive field as a func- 
tion of the wall thickness experimentally observed, with 
a transition between two different modes causing a cusp 
at d^~ T = 13 nm (with d^~ T the thickness at which the 
transition occurs). However, the absolute values com- 
puted for the coercivity are greater than the experimental 
data. 



B. Effect of the stray field 

We ascribe such difference between calculations and 
experimental results to the interaction of each tube with 
the stray fields produced by the array - an effective anti- 
ferromagnetic coupling between neighboring tubes, which 
reduces the coercive field (as previo usly demonstrated in 
the case of nanowires; see Fig. 6a)! 37 " 40 1 In these inter- 
acting systems, the process of magnetization reversal can 
be viewed as the overcoming of a single energy barrier, 
AE. In an array with all the nanotubes initially magne- 
tized in the same direction, the magnetostatic interaction 
between neighboring tubes favours the magnetization re- 
versal of some of them. A reversing field aligned opposite 
to the magnetization direction lowers the energy barrier, 
thereby increasing the probability of switching. The de- 
pendence of the applied field on the energy barrier is 
often described^ by the expression 



AE = U 1 



H 



where H is the applied field, and H n denotes the intrinsic 
coercivity of an isolated wire. For single-domain particles 
having a uniaxial shape anisotropy, the energy barrier 
at zero applied field, U, is just the energy required to 
switch by coherent rotation, K{L). If we assume that 
the switching field H s is equal to H C1 then 



H r = Hi 



Ho, 



(4) 



where H n denotes the intrinsic coercivity H% or H% of 
an isolated tube, and H int corresponds to the stray field 
induced within the array given by 
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first principles, although values between unity and some 
tens could be a reasonable estimate for this quantity.^ 
Besides, E int (s) is the magnetostatic interaction between 
two nanotubes separated by a distance s. Such interac- 
tion can be calculated by considering each tube homoge- 
neously magnetized and is given by 
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The resulting curve H int (d w ) (e = 20) is illustrated in 
the top panel of Fig. 6a. The stray fields produced by an 
array of nanotubes are significant for the experimentally 
investigated tubes, being on the order of 350 Oe for d w = 
5 nm to 580 Oe by d w = 21 nm. 



C. Results 

The hysteresis loop in normalized axis of a sample of 
L = 3 /im, R = 50 nm, and d w = 13 nm is presented in 
Fig. 5. In this loop we have subtracted a paramagnetic 
background. Our results are combined in the lower panel 
of Fig. 6. Experimental data for the coercivity of the 
array are depicted by dots. In this figure we can observe 
the strong dependence of the coercivity as a function of 
the tube wall thickness, evidencing clearly the existence 
of a maximum. Also the coercivities of an isolated tube 
and interacting array obtained from our calculations are 
depicted in the same figure by dashed and solid lines, re- 
spectively. We consider e = 20 in Eq. (4). Note the good 
agreement between experimental datapoints and analyt- 
ical results for interacting arrays for d w > 8 nm. 
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In the previous equation we have assumed that the re- 
versal of individual nanotubes produces a decrease of 
the magnetostatic energy Ei nt that equals the magnetic 
anisotropy barrier AE. Besides, e is an adjustable 
parameter that depends on the distribution of mag- 
netic tubes in space and on the long-distance correlation 
among the tubes. The value of e can not be obtained from 



FIG. 5. (Color online) Magnetic hysteresis loop measured for 
L — 3 /im, R — 50 nm, and d w — 13 nm (after substraction 
of a paramagnetic background and normalization.) 

The deviation of the experimental datapoints from the 
calculated curve for d w < 8 nm likely originates from 
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ferent materials. In the range of parameters considered, 
we observe that an increase of the external radius, R, 
results in an increase of the transition thickness, d^~ T . 
Furthermore, the curves <A^ _T (R) are steeper for mate- 
rials with longer exchange lengths. Figure 7 can also be 
interpreted as a phase diagram, in that each line sepa- 
rates the T mode of magnetization reversal, which pre- 
vails in the upper left region of the (i?, d w ) space, from 
the V mode, found in the lower right area. 
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FIG. 6. (Color online) (Upper panel) Stray field Hi nt calcu- 
lated from Eq. (5) as a function of the tube wall thickness 
d w . As shown in the inset, in an array, contributions from the 
six nearest neighbors, the twelve second-degree neighbors, and 
the tubes situated farthest away from the probe tube combine 
to form the so-called stray field, which is in the — z direction 
and therefore drives the probe tube to reverse its magneti- 
zation at smaller absolute values of the applied field than it 
would if it were isolated. (Lower panel) Coercive field as a 
function of the tube wall thickness for FeaCU nanotubes. The 
red (dark gray) dots correspond to the data measured in or- 
dered arrays, the light blue (dashed) line represents the values 
calculated for an isolated tube, and the deep blue (solid) curve 
is calculated for an array of nanotubes. Parameters: R = 25 
nm, L = 3 /xm and D — 105 nm. Equation (5) was used with 
e = 20. 



the structural imperfections of the tubes. Deposition of 
the magnetic material may lead to granular walls at the 
initial stages of the growth, whereas further increases in 
wall thickness accounted in a smoothening. Other factors 
not account for in our theoretical model include thermal 
instability, which has a stronger effect in thinner parti- 
cles, possible shape irregularities,^ and the finite length 
of the tubes. 

The results presented above may be generalized. We 
now proceed to investigate how the curve will be affected 
by changes in the tube radius, R, and of the material, by 
considering the trajectories of the transition thickness, 
d^~ T . Such trajectories are shown in Fig. 7 for four dif- 



FIG. 7. (Color online) Trajectories of the transition thick- 
ness, dZ~ T , as a function of R, for FesCU (circles, Permalloy 
(squares), Co (triangles) and Fe (stars). Because of the simi- 
larity of Ni and FesCU in terms of their magnetic parameters, 
the results presented for FesCU approximate the case of Ni 
tubes as well. 



IV. CONCLUSIONS 

In conclusions, by means of simple models for the do- 
main boundary that appears during the magnetization 
reversal in nanotubes, we can calculate the coercive field 
in ordered arrays of ferromagnetic nanotubes as a func- 
tion of the tube wall thickness and the radius. A transi- 
tion between two different modes of magnetization rever- 
sal, from a vortex boundary, in thin tubes, to a trans- 
verse boundary, in thick tubes, is responsible for the 
non-monotonic behavior of the coercivity as a function 
of wall thickness experimentally observed. The effect of 
the stray field originating from the magnetostatic inter- 
actions between the tubes of the array must be included 
to obtain a quantitative agreement between experimen- 
tal and theoretical results. Because of its long range, 
the magnetostatic interaction strongly influences the co- 
ercivity of the array. Finally, the presence of a coercivity 
maximum at a certain optimum wall thickness should be 
a quite general phenomenon, observable for a variety of 
ferromagnetic materials and of tube radii. Experimen- 
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tal work remains to be done in order to validate these 
predictions. 
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